arXiv:1509.01422v3 [cond-mat.mtrl-sci] 17 Nov 2015 


Theory of Curie temperature enhancement in electron-doped EuO 

Tobias Stollenwerk^’^ and Johann Kroha^’^’B 

^ Physikalisches Institut and Bethe Center for Theoretical Physics, 

Universitdt Bonn, Nussallee 12, 53115 Bonn, Germany 
^German Aerospace Center, Linder Hohe, 51147 Cologne, Germany 
^Center for Correlated Matter, Zhejiang University, Hangzhou, Zhejiang 310058, China 
(Received 10 September 2015; revised manuscript received 20 October 2015; published 17 November 2015) 

We present a comparative, theoretical study of the doping dependence of the critical temperature 
Tc of the ferromagnetic insulator-metal transitions in Gd-doped and 0-dehcient EuO, respectively. 
The strong Tc enhancement in Eui_a;GdxO is due to Kondo-like spin fluctuations on the Gd sites, 
which are absent in EuOi-a;. Moreover, we hnd that the Tc saturation in Eui-^Gd^O for large 
X is due to a reduced activation of dopant electrons into the conduction band, in agreement with 
experiments, rather than antiferromagnetic long-range contributions of the RKKY interaction. The 
results shed light on possibilities for further increasing Tc- 

PACS numbers: 71.30.-|-h, 75.20.Hr, 75.30.—m, 75.50.Pp 


I. INTRODUCTION 

The demand for ever increasing speed and integrabil- 
ity of magnetic information storage devices as well as 
other spintronics applications calls for materials that are 
capable of transforming electronic or optical signals ef¬ 
ficiently into magnetization and vice versa. Electron- 
doped europium monoxide (EuO) is a promising candi¬ 
date for this purpose, as it undergoes a simultaneous fer¬ 
romagnetic (EM) and insulator-to-semimetal transition 
[l|, exhibiting an outstanding magneto-electric response, 
includiM the strongest colossal magnetoresistance effect 
known , magneto-optical effects , and a spin po¬ 
larization of the conduction band of nearly 100 % in the 
EM state d0|- Improved sample fabrication techniques 
mm and europium monoxide’s epitaxial integrability 
into Si 0 and GaAs [l^ structures have renewed and 
intensified the interest in this material during the past 
few years. 

Stoichiometric EuO is a wide band gap semiconductor 
with rocksalt structure. The local magnetic moments of 
m = 7/2 11B situated in the Eu 4/ shell constitute a pro¬ 
totype Heisenberg ferromagnet, ordering ferromagneti- 
cally at the Curie temperature of Tc = 69 K. Their inter¬ 
action is mediated by virtual excitations (hybridization) 
of the tightly bound Eu 4/ electrons into the spatially 
more extended Eu 5d orbitals and an exchange between 
the latter [l^. Upon electron doping the EM transition 
is accompanied by a simultaneous insulator-to-semimetal 
transition with a resistivity drop by 8 to 13 orders of 
magnitude [iilil. Raising the transition temperature 
significantly is one of the major challenges involved in 
bringing the extraordinary properties of doped EuO to 
technological use. 

Since early on, the general trend of Tc being increased 
by electron-doping has been associated with the forma¬ 
tion of magnetic polarons (l^fT^ . i.e., conduction elec¬ 
trons dragging along a magnetic polarization cloud of 
local Eu 4/ moments. However, the experiments re¬ 
veal more complex behavior. Gadolinium doping re¬ 


places Eu by Gd atoms, leaving the lattice of magnetic 
4/ moments intact and donating one additional elec¬ 
tron per Gd atom from the Gd 5d shell. In the doped 
material, Eui_a;Gda;0, Tc increases to values between 
120 K and 170 K for Gd concentrations in the range 
of a; = 0.04,..., 0.08 (^ fisl - l^ . depending on the sam¬ 
ple quality, strain, and measurement conditions (^ . In¬ 
variably, Tc saturates for higher x. Oxygen defects in¬ 
troduce two electrons per O defect, but only a weak 
7T increase has been observed in bulk, Eu-rich EuOi-a, 
. The Tc increase reported in Ref. [2^ 
for EuOi-x may presumably be attributed [22| to the 
presence of a large external magnetic field inherently nec¬ 
essary for the SQUID measurement technique used. 

The magnetic polaron theory alone cannot account 
for the Tc saturation at high Gd concentration nor for 
the fact that O defects essentially do not raise Tc, even 
though they introduce twice as many carriers per impu¬ 
rity as Gd doping. It has been proposed that the 

Tc saturation might be understood in that, for increasing 
conduction-band filling, the oscillatory RKKY interac¬ 
tion [28l - [^ acquires increasingly antiferromagnetic (AE) 
contributions due to the decreasing RKKY wavelength. 
This requires, however, unrealistically high band filling. 
On the other hand. Hall resistivity measurements indi¬ 
cate ( 2 ^ that the density of mobile charge carriers acti¬ 
vated into the conduction band saturates in line with the 
Tc saturation, providing a phenomenological reason for 
the limited Tc increase in Eui_a;Gda;0. In theoretical cal¬ 
culations, treating the O vacancies in a static approxima¬ 
tion, Sinjukow and Nolting found no increase of Tc in 
EuOi-a; for an appropriate choice of system parameters. 
More sophisticated resummation techniques for the Gd 
impurity potential in Eui_a;Gda;0 were able to correctly 
describe a shallow maximum of Tc as a function of Gd 
concentration (^ . but not the saturation of the mobile 
charge carrier density Recent ab initio calculations 
provided more quantitative results on the coupling 
parameters and spectral densities, but did not lead to 
a consistent understanding of all the experimental facts 
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described above. Taking the local spin fluctuations on 
Gd impurities into account, Arnold and Kroha could 
explain details of the magnetization behavior, the simul¬ 
taneity of the insulator-seniimetal transition and the Tc 
increase in Eui_a;Gda;0. 

In this paper we report a comprehensive, theoretical 
study of the Tc enhancement in electron-doped EuO, ex¬ 
tending the model of Ref. . Gd impurities as well as 
O defects are treated dynamically as Anderson impuri¬ 
ties hybridizing with the conduction band, however with 
strong or moderate on-site repulsion, respectively, ensur¬ 
ing the single or double occupancy of the Gd impurity 
or O vacancy orbitals. The direct comparison of the two 
cases shows that indeed the Tc increase with Gd dop¬ 
ing is caused by the Kondo-like spin fluctuations on the 
Gd sites and the concatenated accumulation of impurity 
spectral weight as well as conduction electron spectral 
weight at the chemical potential. This dynamical many- 
body effect drives the metallic transition, which, in turn, 
enhances the polaronic EM coupling between the Eu mo¬ 
ments. By contrast, on O vacancies the two defect elec¬ 
trons form a spin-singlet, and local spin fluctuations are 
absent, leading only to a moderate Tc enhancement due 
to a weak population of the conduction band. Moreover, 
the theory explains that in Eui_xGda;0 the activation of 
charge carriers into the conduction band decreases with 
increasing doping concentration, in agreement with ex¬ 
periments 1^, leading to the saturation of Tc- Includ¬ 
ing explicitly the RKKY interaction in our theory, we find 
that for all relevant temperatures T and doping concen¬ 
trations (band fillings) its long-distance AF contributions 
are negligible. 

The paper is organized as follows. In Sec. II, we give 
a detailed justification of our model for Eui_xGda:0 and 
EuOi-x and describe the theory for its evaluation. The 
results are shown and discussed in Sec. III. We conclude 
in Sec. IV with a suggestion for a possible pathway to 
further enhance the transition temperature in electron- 
doped EuO. 


II. THEORY 
A. Model 

The model Hamiltonian for Eui-a^Gd^O as well as 
EuOi-x reads, 

H = HQ + H,f + H,mp ■ ( 1 ) 

The conduction band, comprised of the hybridizing Eu 
5d6s orbitals, is described by 

Ho = ^(ek - M)cLck^ > (2) 

kfT 

where is the conduction electron creation operator 
and Ek the conduction band dispersion, /i is the chemical 
potential that fixes the total electron number (conduction 


and impurity electrons). In undoped EuO, p lies in the 
gap below the conduction band. The lattice of Eu 4/ 
moments is described by a Heisenberg Hamiltonian and 
a local coupling term between the Eu 4/ moments and 
the conduction electron spins, 

Hcf = -J2 ■ S, - Jcf ■ (3) 

d.i) * 

Here, is the 4/ spin, ms = —7/2,...,-1-7/2, and 

(Ti = i X^CTcr' ''L'^’o-o-'Ocr' couduction electron spin 

operator at site i. Jij > 0 is the direct exchange coupling 
between the localized moments which is independent of 
the conduction band occupation and therefore responsi¬ 
ble for the Gurie temperature of 69 K in stoichiometric 
EuO. Jcf is the exchange coupling between the 4/ mo¬ 
ment and the conduction electron spin ai at lattice site i. 
The Gd impurities and O vacancies are described as An¬ 
derson impurities with a single electron binding energy 
Ed < Q and a hybridization with the conduction band, 

Hcd = {Ed-fi)J2dld^a (4) 

fil'T {j} 


where is the electron creation operator in an impurity 
or defect orbital at site j and {/} indicates a summation 
over the randomly placed impurity sites. In the following 
we will use the term impurity for both, Gd impurities 
and O vacancies, unless stated otherwise. The impurity 
number density will be denoted by nj. Gd carries one 
extra electron in the hd shell as compared to Eu. Hence, 
Gd is in the strongly correlated regime with a strong 
onsite repulsion U preventing double occupancy. Due to 
stoichiometry, an O vacancy attracts two extra electrons 
from the surrounding metal ions. Therefore, it is in the 
weakly correlated regime, with double occupancy in the 
ground state, i.e., a moderate onsite repulsion 0 < 17 <S 
\Ed\. The model, Eqs. ffl-ffl, inherently incorporates 
the RKKY interaction |28l - l30l | via 2nd-order, non-local 
perturbation theory in Jcf. Since, apart from RKKY 
effects, only local self-energies will be important for the 
following treatment of this paper, we here include the 
RKKY Hamiltonian explicitly. 


H 


RKKY_ 


E 

(iJj) 


+ (SfS^ + SfSy)] .(5) 


It is to be amended to the model Hamiltonian O- A rec¬ 
ollection of its derivation and the expressions for the cou¬ 
pling constants are given in Appendixsee Eqs. (IA8II . 
(lAl) . 


B. Selfconsistent theory 

To evaluate this model, we follow and extend Ref. . 
While the large spins of the 4/ Heisenberg lattice can 
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FIG. 1: (Color online) Conduction band (upper panel) and 
Gd impurity (lower panel) spectral density across the phase 
transition in Eui_a;Gda;0 for x = 0.01. Tc « 95 K. Solid 
curves represent the majority (a ='() and dashed curves the 
minority (a =4-) spin spectral density. 


FIG. 2: (Color online) Conduction band (upper panel) and O 
vacancy (lower panel) spectral density across the phase tran¬ 
sition in EuOi-a; for x = 0.01. The bare O-defect parameters 
are Ed ~ 0.0Do, U = 0 (double occupancy). Tc ~ 78 K. 
Solid curves represent the majority (a =t) and dashed curves 
the minority (a =4') spin spectral density. 


be treated in mean-field theory, it is essential to de¬ 
scribe the Anderson impurities dynamically, in order to 
account for their spin and charge fluctuations. The con¬ 
duction electron selfenergy induced by the impurities will 
be treated in a single-site approximation, i.e. it is given 
by the full impurity T matrix times the impurity con¬ 
centration m. This is valid for dilute impurities, where 
inter-impurity correlations are negligible. Writing the 
(retarded, uj = uj + iO) conduction electron Green’s func¬ 
tion as, 


GcaiK^^) 


1 

U} + ^ - Sk - Sccr(a;) 


( 6 ) 


states (DOS). Jif = •^oj is the effective mean-field 
coupling of the Heisenberg lattice. The short-range 
nature of Jij restricts the lattice summation essentially 
to the nearest neighbors of site z = 0. The mag¬ 
netic polaron effect (l^H^ is incorporated in Eq. ([5]) 
via the conduction electron magnetization {a). For 
later use, the conduction carrier density is given by, 
ric = / dujf{uj)Aca{u!)- The calculation of the local 

impurity Green’s function Gda-i<-^) depends on whether 
the impurity is in the strongly (Gd impurities) or the 
weakly (O vacancies) correlated regime. 


yields for the total conduction selfenergy, 

^ca{^) =niV‘^Gda{^) - crJcfiS), CT = ±i, (7) 

where Gda{'^) is the full Green’s function of the impurity 
electrons. The expectation value of the 4/ spins is de¬ 
termined in mean-field theory by (with /3 = l/ksT, the 
inverse temperature), 

7 5'e'®(2J4/(S>+J,/<<7))s 

(*5) = -, (8) 

y-2^ ^ gl3{2Jif{S) + J,f{a))S 

Here, the conduction electron magnetization (cr) is cal¬ 
culated as, 

(cr) = ^ y dufiuj) [Act(w) - Ac^(a;)] , (9) 

with Accr(uj) = — i ImGccr(k, w -b zO) the spin- 
dependent, interacting conduction electron density of 


a. Gd impurities. We employ the auxiliary particle 
technique in non-crossing approximation |35l - j37| to de¬ 
scribe the spin and charge fluctuations in the Gd 5d or¬ 
bitals. The limit f/ —oo is taken, for simplicity, in order 
to prevent double occupancy. Since in Eui_2,Gda;0 the 
DOS near the chemical potential is low or even vanishing, 
the Hondo temperature of the fluctuating spin, Tj^, is far 
below Tc- The NGA is known to give reliable results for 
energies above and down to somewhat below Tr-. In a 
magnetic field it produces, in addition to the Zeeman- 
split Hondo resonance, a spurious low-temperature sin¬ 
gularity at the Fermi level for T < Tk- However, since 
Tff <C Tc in Eui_a,Gda;0, the effect of this singularity is 
negligible for the temperature range relevant here. The 
NGA is also versatile enough to include an arbitrary en¬ 
ergy dependence of the DOS. Therefore, it is the appro¬ 
priate method for the present purpose [sj. The NGA 
equations, adapted for the Eui-^Gd^O system with a 
gapped DOS and a non-trivial chemical potential, are 
given in Appendix iBl 
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b. 0 vacancies. The weak interaction on the O de¬ 
fects, where spin fluctuations are negligible, can be ac¬ 
counted for in the second-order perturbation theory in 
U. The (retarded) 0-defect electron Green’s function is 
Gdo-(w) = 1/ [uj + — Ed — Sdcr(w)], and the correspond¬ 
ing selfenergy reads, Sdo-(w) = (uj) + (w), where 

= 7tV^A,^{uj) + U J def{e)Ad-Ae) (10) 

ImE("j(cc) = (11) 

-U'^ J dei J de2 Ada{ei + oj)Ad-a{ei + e2)Ad-aie2) 

X IKsi) + /(^i + w)] [/(e 2 ) - f{ei + 62)] , 
where /(w) and b{u;) are Fermi and Bose distribution 

/q\ 

functions, respectively, and (uj) is given by the 

Kramers-Kronig relation. 

Long-range RKKY coupling. When the RKKY inter¬ 
action is included, the 4/ Heisenberg coupling is changed 
to Jij —> Jij -|- Kij . In the above equations this leads to 
a modified mean-field coupling, 

•/4/= E (12) 

jn.n.Q 

The lattice summation in the second term is carried out 
over the fee lattice of the EuO rocksalt structure. Note 
that on mean-field level only the longitudinal RKKY 
component contributes and can give FM as well as AF 
contributions to the total coupling. involves non¬ 

local Green’s functions (c.f. Appendixl^ and, hence, the 
band dispersion et- For simplicity and since anisotropy 
effects are not important in bulk EuO, we assume for the 
RKKY interaction an isotropic dispersion which is con¬ 
structed such [ 3 ^ that it reproduces the bare conduction 
DOS. 



FIG. 3: (Color online) Curie temperature vs. the doping con¬ 
centration m for Gd-doped and O-deficient EuO. Bare pa¬ 
rameters for O vacancies: Ed = 0.0Do, U = 0. The blue 
curve represents the behavior for Eui_ 2 :Gda ;0 including the 
RKKY interaction, see text. 



FIG. 4: (Color online) Curie temperature vs. on-site Coulomb 
repulsion U in EuOo.99 for various impurity level energies Ed. 
It is seen that for more tightly bound defect electrons the Tc 
enhancement in EuOa-x is even weaker than for the param¬ 
eter values of Fig. |3] 

The system is subject to the doping condition that 
the total density of charge carriers is ntot = n-i for 
Eui_a;Gda;0 and ntot =2ni for EuOx- That is, 

E / dujf{uj) [Acct(w) -b njAdcriuj)] - ntot = 0 . (13) 

<T d 

The selfconsistent set of equations (l^- (fT51) . in the case 
of Eui_a;Gda;0 amended by the NCA equations (IBlI) - 
(IB3I) . is solved by iteration, where in each iteration step 
the chemical potential p is adjusted so as to fulfill the 
particle number constraint (ITHl) . Note that the RKKY 
coupling strength , Eq. (IA9I) . is also determined self- 
consistently via the interacting conduction electron prop¬ 
agators. 


III. RESULTS 
A. Parameter values 

For the numerical calculations below, we 
choose a semi-elliptic DOS for the non-interacting 
conduction band of stoichiometric EuO, {uj) = 
( 2 /^Do) — ptQ — DqY/Dq — The conduction 

half-bandwidth is taken to be Dq = 8 eV, and the chemi¬ 
cal potential of the undoped system lies in the gap below 
the conduction band, /ig = — 0 . 02 Do, consistent with ex¬ 
periments 0. All energies are measured relative to the 
(interacting) chemical potential p and are given in units 
of Dq. The mean-field Heisenberg coupling J 4 y (without 
RKKY interaction) is chosen such that Tq = 69 K is ob¬ 
tained for undoped EuO. This yields J 4 / = 7 • 10“^ Dg 
( 3 ^ . Jof is much larger than J 4 ^, because the overlap 
between the neighboring Eu 4/ orbitals is much smaller 
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FIG. 5: (Color online) Longitudinal RKKY coupling /S'll(r) 
for X = m — 0.04 and various temperatures in Eui_a;Gda;0. 
The inset expands the AF region around two to four lattice 
spacings a of the FCC lattice. 


than their overlap with the conduction wave functions. 
From the spatial separation of the Eu 4/ orbitals the ratio 
Jcf/J^f is roughly estimated to give Jc/ = 0.05 £>0, see 
also Ref. [sj. This also determines the RKKY coupling 
strength via Eq. (IA9I) . We fix the bare parameters of a 
Gd impurity such that for T = 0 and vanishing impurity 
concentration its occupation is ~ 1 and that the impu¬ 
rity electron gets thermally activated in the experimen¬ 
tally relevant temperature range. This yields, Ed = 0.0, 
r := ttV'^IDq = 0.05 £> 0 , and U ^ oo. Note that hy¬ 
bridization and interaction substantially renormalize the 
impurity level, Ed ^ Ed ^ —0.02Do (Haldane shift (3^ . 
see also Fig. [1]), so that ~ 1 is realized in the n/ —>■ 0 
limit. This also renders the Tk of the impurity far below 
Tc, since in our system the DOS at the Fermi level Ep 
remains always Aco-(O) <C 1/Do- For O vacancies, in ab¬ 
sence of more detailed information about their structure 
other than double occupancy, we set the effective hy¬ 
bridization r = 0.05 Do, as for Gd, and perform a scan 
oi Ed < 0 and U within the bound-state, weakly corre¬ 
lated regime, see below. The RKKY interaction will be 
included and discussed in Subsec. IlII Cl only. 


B. Gadolinium impurities vs. Oxygen vacancies 

Figuret[T] and [2] show the evolution of the conduction 
band and impurity spectral densities across the phase 
transition for low-doped Eui_ 2 ,Gd 3,0 and EuOi-x, re¬ 
spectively. For both, Eui-^Gd^O and EuOi-x, in the 
high-temperature insulating phase the spin degenerate 
conduction DOS is comprised of a large, unoccupied band 
and a small side band which is induced by the hybridiza¬ 
tion with the impurity orbitals and is centered around 
the impurity binding energy Ed, lying entirely below /i 
and, therefore, not contributing to the conductivity. As 
the temperature is lowered, in Eui_xGda;0 the impu¬ 



FIG. 6: (Color online) Total magnetization M = (S) -I- 
(a) vs. temperature T for various doping concentrations in 
Eui_a;Gda,0 with long range RKKY coupling (solid curves) 
and without RKKY coupling (dashed curves). 


rity spectrum accumulates spectral weight at the chem¬ 
ical potential which eventually develops into a peak at 
w = 0 (Fig. [1]). Below Tc the spectral densities are split 
into majority and minority bands. Our NGA calculations 
show that this is due to local, Kondo-like spin fluctua¬ 
tions in the Gd hd orbitals [s^ . Because of hybridization, 
the conduction electron DOS develops spectral weight at 
w = 0 as well, and the side band merges with the main 
conduction band. This drives the metallic transition and 
simultaneously enhances, via the magnetic polaron effect 
[c.f. Eq. ([5])], the FM transition temperature as well. 
In EuOi-a;, the local spin fluctuation effect is absent. 
Here, the metallic transition occurs only when the con¬ 
duction side band is eventually broadened and shifted, 
via hybridization with the O vacancy band enough to 
gain overlap with the chemical potential (Fig. [21), lead¬ 
ing to a much lower Tc than in Eui_a;Gda:0. In Fig. [3] 
the doping-dependent Tc enhancement is compared for 
magnetic Gd impurities and non-magnetic O vacancies 
(black and red curves). Here, for O vacancies, U = 0 
(double occupancy) and otherwise the same parameter 
values as for Eui_2;Gda;0 are assumed. This allows for 
a direct assessment of the importance of low-lying, lo¬ 
cal spin fluctuations for the Tc enhancement. The es¬ 
sential role of on-site correlations as well as conduction 
electron doping is further substantiated by Fig. |4l where 
Tc is shown for varying U and Ed values in EuOi-a,: Tc 
is enhanced by repulsive onsite correlations (increasing 
U) and is reduced by the dopant electrons more tightly 
bound to the defect (more negative Ed). 


C. RKKY interaction in Eui-iGdj,© 

We now study the influence of the long-range RKKY 
interaction on the phase transition in Eui-^Gd^O. For 
short Fermi wavelength, the RKKY interaction might 
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make an AF contribution to the total coupling and, thus, 
lead to the experimentally observed saturation of Tq 
(^ flsl - E^ . as has been suggested in Ref. [ 2 ^. In order to 
analyze the possible influence of the RKKY interaction 
on the saturation at high doping concentration, we now 
adjust the value of Jcf such that the theory including 
RKKY reproduces the previous results without RKKY 
interaction (Sec. IlII Bl) in the low-doping regime, and 
will compare the results at high doping. This yields the 
new value J^f = 0.0405 Z?o- The RKKY coupling K^\r), 
selfconsistently calculated for the interacting system, is 
displayed in Fig.[S]as a function of distance r for a typical 
Gd doping concentration of n/ = 0.04 over the complete, 
relevant T range. It shows weak AF behavior only in 
the range of about 2 to 4 FCC lattice constants. The re¬ 
sulting total magnetization M is shown in Fig. [51 While 
the FM magnetization is even enhanced by AT 11 below the 
transition, it does not substantially alter Tc, especially 
for higher doping. This is plausible, because the RKKY 
interaction is not active for T > Tc (empty conduction 
band), but its long-range, overall FM behavior enhances 
M once the band is populated for T < Tc ■ Such enhance¬ 
ment of the FM coupling by an RKKY-like interaction is 
consistent with recent experiments on EuO doped with 
non-magnetic La atoms [^. Fig.|2]directly compares Tc 
with and without RKKY coupling (blue and black curves, 
respectively) in our calculation. It is seen that including 
the RKKY interaction and reducing the direct exchange 
coupling to Jcf = 0.0405 Dq not only reproduces the Tc 
behavior at small nj (by construction), but also does 
not change the behavior for the largest n/ considered. 
For the small band fillings relevant in Eui-^Gd^O the 
effects of the RKKY interaction are essentially doping 
independent and can be absorbed in a proper choice of 
Jcf, at least as far as Tc in bulk systems is concerned. 
We conclude that the experimentally observed Tc sat¬ 
uration behavior in Eui_a;GdxO for large n/ cannot be 
attributed to the RKKY interaction. Note that ab ini¬ 
tio calculations Haul presumably overestimate the anti¬ 
ferromagnetic contributions from the RKKY interaction, 
because they do not take the Kondo-like spin fluctuations 
on the Gd sites and the resulting accumulation of spec¬ 
tral weight at the chemical potential into account. As a 
consequence, the RKKY wavelength comes out too short 
and, hence, its antiferromagnetic contributions too large. 
This may be the origin why these calculations overesti¬ 
mate the decrease of Tc for large doping concentration 
as compared to experiments pm. 


D. Dopant activation and Tc saturation 

In Fig. 13 we show the charge carrier concentration in 
the conduction band Uc (number of carriers per lattice 
unit cell; inset) as well as the dopant activation Uc/ni as 
a function of impurity concentration nj at the lowest tem¬ 
perature considered, T = 5 K, similar to Ref. For 
low doping the impurity spectral weight Adai^) (both, 



FIG. 7: (Color online) Dopant activation Uc/nj for Gd doped 
and O deficient EuO at the lowest temperature considered, 
T = 5K. The carrier concentration Uc is shown in the inset. 


cr =t, i) is almost entirely shifted above Ep at this tem¬ 
perature, as can bee seen from Figs. (TJ |3 lower panels, 
so that the impurity level is completely emptied into the 
conduction band. Consequently, the dopant activation is 
nc/nj = 1 (Gd) or Uc/nj = 2 (O vacancies) up to a dop¬ 
ing concentration of about nj = 0.01 (Fig.|3). For higher 
71/, the impurity contribution to the conduction electron 
selfenergy, Eq. gets increasingly enhanced by the dis¬ 
order. Via the Kramers-Kronig relation for the real part 
of I]cct(w) this implies a downward shift of the conduc¬ 
tion side band and, connected with it, a downward shift 
of the impurity band below Ep. This is seen in Fig. |8l 
It results in a re-population of the impurity levels and a 
crossover to a reduced ric/n/, as seen in Fig.|3 Note that 
the description of the Gd impurities as Anderson impu¬ 
rities with spin fluctuations is crucial for the downward 
shift of the impurity levels. The reduction of the dopant 
activation Uc/ni is in agreement with the experimental 
findings of Ref. 0- Note that in Ref. a a reduced 
dopant activation is also found in the limit of small np 
Presumably this is, because their Hall measurements de¬ 
termine the mobile carrier density n only. However, for 
small impurity concentration, part of the electrons in the 
conduction band will be bound around the impurity lo¬ 
cations. However, all electrons in the conduction band, 
given by ric, contribute to the electron-enhanced mag¬ 
netic coupling. Gomparing the doping dependence of Hc 
in the inset of Fig. |3 with Tc in Fig. |3| shows that the 
latter follows the behavior of Uc- Displaying now Tc 
(same data as in Fig. 131) as a function of the carrier con¬ 
centration Uc in Fig. [9] shows that it grows essentially 
linearly with Uc, showing only a slightly decreasing slope 
for the highest Uc- Note that the highest carrier concen¬ 
tration of Uc ~ 0.04 appears experimentally achievable, 
while the corresponding doping concentration of nj = 0.4 
is not, due to stability reasons of the crystal structure. 
The saturation-like behavior of Tc for large nj doping 
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FIG. 8: (Color online) Majority spectral densities for conduc¬ 
tion electrons (upper panel) and impurity electrons (lower 
panel) in Eui_a;Gda;0 at T = 5 K for various impurity con¬ 
centrations ni. The figure shows the downward shift of the 
impurity levels with increasing nj. 


is, therefore, to be considered a consequence of the re¬ 
duced dopant activation for large do ping , in complete 
agreement with the conclusion of Ref. . 



FIG. 9: (Color online) Curie temperature vs. carrier con¬ 
centration ric in the conduction band for Gd-doped and O- 
deficient EuO. (Inset) Semilogarithmic plot of the data. 

generating low-energy spin fluctuations and the charge¬ 
doping impurities need not necessarily be of the same 
type. More efficient carrier doping may be achievable by 
using different types of donor atoms (with impurity levels 
closer or above the Fermi energy) in addition to Gd, or 
by carrier coupling at interfaces in heterostructures. 


IV. CONCLUSION 

We have performed a detailed comparison of the 
FM insulator-metal transitions in Eui-^Gd^O and in 
EuOi-x, respectively, using a model that treats the 
dopant impurities as Anderson impurities in the strongly 
(Gd) or weakly (O vacancies) correlated regime, and that 
had previously provided a detailed description of ex¬ 
perimental magnetization, resistivity and total conduc¬ 
tion band polarization data. Our results show that for a 
significant, doping-induced Tq enhancement a coopera¬ 
tion of two effects is necessary, (1) Kondo-like, low-energy 
spin fluctuations accumulating impurity as well as con¬ 
duction spectral weight at the Fermi energy and (2) ef¬ 
ficient population of this low-lying spectral weight and 
subsequent enhancement of the FM interaction between 
the 4/ moments mediated by the conduction electrons. 
In addition, our calculations provide evidence that the 
tendency of Tq to saturate for high doping concentra¬ 
tions is not due to AF contributions of the RKKY in¬ 
teraction but rather due to a limitation of the dopant 
electron activation into the conduction band, confirming 
experimental results Hence, an increase of Tq be¬ 

yond the presently achievable values may be possible, if 
only the conduction band can be populated in a more 
efficient ww. This is in line with recent pump-probe ex¬ 
periments [3 where enhanced FM coupling was achieved 
by photodoping into the conduction band |4l| . The com¬ 
bination of all these findings point to a possible pathway 
to further enhancement of Tc'. the magnetic impurities 
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Appendix A: RKKY interaction in paramagnetic 
metals and semimetals 

In this appendix we derive the expressions for the 
RKKY interaction, allowing for an arbitrary magneti¬ 
zation of the conduction band and then specializing for 
the paramagnetic case (vanishing magnetization) and the 
semimetallic case (complete magnetization). The inter¬ 
action Hamiltonian between localized spins at sites i 
and the conduction electron spins reads, 

Hcf ~ Jcf ^ ^ Sz ' , (-^1) 

i 

where = 1/2 conduction electron 

spin operator at site i and cr = (cr^, cr^) the vector 
of Pauli matrices. Evaluating the time evolution of the 
conduction electrons in the presence of another localized 
spin Sj according to exp[— dTHcfir)] in first order 
of the spin coupling Jcf and tracing out the conduction 
electron degrees of freedom, one obtains in the static limit 



























(w = 0) the RKKY interaction operator between the local 
spins Si and S^, 

= -J^cf ((S. • s,)(S, • . (A2) 

Here ((... ))^ := trc{e“^^(... )}/Zc denotes the thermal 
trace over the conduction electron Hilbert space. Using 
Wick’s theorem, it can be written as, 


tjRKKY 



OL,0—x,y,z (Try' 


(A3) 


where H^?’ is the conduction electron density propagator 
between the sites i and j as depicted diagrammatically 
in Fig. [ini It has the general form, 

^ij ~ '^ GjityjiCn ~\~ iuj')Gij a'ij'^n) • (A4) 



FIG. 10: Diagram for the spin-dependent conduction elec¬ 
tron polarization function 11^°^ (a;), generating the RKKY in¬ 
teraction. The solid lines represent conduction electron prop¬ 
agators. 

interaction is in general anisotropic for a magnetized con¬ 
duction band. 

We now present explicitly the expressions for the spe¬ 
cial cases of a paramagnet and of a semimetal. For a 
paramagnetic conduction band we have Gij a = Gij^-a 
independent of spin. Hence, the RKKY coupling is 
isotropic, and we have the paramagnetic RKKY Hamil¬ 
tonian, 


In the static limit it reads, 

= - J de fie) X (A5) 

[^ijcTi£)^GGij a' (e) Aij^ii£)^eGij„ief\ , 

where Aij ^{e) = —ImG^-cr(e + *0)/7r. Performing the 
spin contractions in Eq. (|A3I) and defining the longitu¬ 
dinal and the transverse polarization functions, respec¬ 
tively, as 

ny(0) = ^En?,"(o) (A6) 

a 

= “E de fie) Aij„i£)ReGijaie) 

a 

nq(o) = ^Enr"(o) (A7) 

(7 

= “E / de fie) Aij „ie)ReGij -„ie) , 

one obtains the RKKY interaction Hamiltonian, 

H^kky ^ (A8) 

= - E K + s^s^) 

(i,j) 

where the sum runs over all (arbitrarily distant) pairs of 
localized spins and Sj, and 

ipii = (0), Kt, = (0), (A9) 

are the longitudinal and transverse RKKY couplings, re¬ 
spectively. As seen from Eqs. (IA8I) and (IA9|I . the RKKY 


= (AlO) 

(hi) 

with 

= -^E / de fie) A,j „ie)ReG^j „ie) . (All) 

For a semimetal, i.e., for a completely spin-magnetized 
conduction band with majority spin a =t we have 
Aij|(e) = 0, and the semimetallic RKKY Hamiltonian 
reads, 

= (A12) 

(hi) 

^ {sfs^^+sjsy)] , 

with 

= -^/ defie)A,,^ie)ReG,,^ie) (A13) 

^ ^ j de fie) A,j^ie)ReG^j ^ie) . (A14) 

The missing spin summation in Eqs. (jA13l) . (IA14I) as 
compared to Eq. (lAlip indicates that in the completely 
magnetized band only the majority spin species con¬ 
tributes to the coupling. Note, however, that the trans¬ 
verse coupling JpMY^ij i® non-zero even in the ferro- 
magnetically saturated case because of virtual (off-shell) 
minority spin contributions represented by the real part, 
ReGy 4 ,(£) in Eq. (|A14I) . The Curie temperature Tc, 
where the band magnetization vanishes, is controlled by 
the paramagnetic RKKY coupling, Eq. (lAlOl) . while far 
below Tc the carrier concentration nc is determined by 
the semimagnetic RKKY coupling, Eq. (jA12l) . 
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Appendix B: NCA equations 

The selfconsistent NCA equations for the pseud¬ 
ofermion (/) and slave boson (6) self-energies read 
[uj = uj -\- iO) 

Ef^{uj) = V^ J de[l - f{e)]A,^ie)Gb{uj - e) (Bl) 

= f dsfie)Aca{e)Gf4uj + e), (B2) 

with the auxiliary particle Green’s functions, 
Gfa{uj) = 1/[ui + \-Ed-T.fa{uj)] and 

Ghioj) = 1/[w — A — Eb(a;)], respectively. A is a 

positive parameter, taken to A —^ cxd in order to effect 


the constraint on the auxiliary particle number operator, 
Scr= 1- Note that these NCA equations 
are coupled to the equations ([H])-© for the interacting 
conduction electrons via the common chemical potential 
p and via the conduction electron DOS of the interacting 
system in presence of a dilute, but finite impurity con¬ 
centration, Aco-(e)- The Gd impurity electron Green’s 
function is obtained from G/o-, Gb as, 

Gda{^) = J -^[Ab{£)Gfa{e + uj) - Afa{e)Gl{e - uj)] . 

(B3) 

For an efficient and accurate method for numerically solv¬ 
ing the set of equations (IBl() - (jB3l) see Ref. [i^ . 
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